Spin relaxation near the metal-insulator transition: dominance of the Dresselhaus spin-orbit 

coupling 
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We identify the Dresselhaus spin-orbit coupling as the source of the dominant spin-relaxation mechanism 
in the impurity band of doped semiconductors. The Dresselhaus-type (i.e. allowed by bulk-inversion asym- 
metry) hopping terms are derived and incorporated into a tight-binding model of impurity sites, and they are 
shown to unexpectedly dominate the spin relaxation, leading to spin-relaxation times in good agreement with 
experimental values. This conclusion is drawn from two complementary approaches employed to extract the 
spin-relaxation time from the effective Hamiltonian: an analytical diffusive-evolution calculation and a numeri- 
cal finite-size scaling. 

PACS numbers: 72.25.Rb, 76.30.Pk, 72.20.Ee, 03.67.-a 
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Spin dynamics in semiconductors is a fundamental issue in 
view of the rich physics involved and the potential techno- 
logical applications lUlQl- It is thus not surprising that spin 
relaxation studies were already performed in the early days of 
semiconductor research and are intensely pursued today 
with modern experimental techniques fllll- An intriguing ex- 
perimental observation is the fact that in n-doped semiconduc- 
tors at low temperatures the spin relaxation time Ts presents a 
maximum as a function of the doping density near the metal- 
insulator transition (MIT) ill lei IsUlJO . 

Interestingly, while the mechanisms behind spin relaxation 
have been properly identified at high ternperatures or for dop- 
ing densities away from the critical one a theoret- 
ical understanding of low-temperature spin relaxation close to 
the MIT is still lacking. This unsatisfactory state of affairs has 
motivated some attempts to identify the relevant mechanisms 
for spin relaxation [,l5r; l81 close to the MIT. In particular, on 
the metallic side of the transition, Shklovskii has proposed the 
applicability of the well-known Dyakonov-Perel mechanism 
usually valid in the conduction band |15]. Furthermore, a 
tight-binding model of impurities including Rashba spin-orbit 
coupling has been developed |18]. The spin relaxation times 
resulting from this last model were larger than the experimen- 
tal values, implying that other mechanisms should be active in 
this density regime. 

In this work we identify the Dresselhaus coupling as the 
source of the leading spin relaxation mechanism on the metal- 
lic side of the transition. This conclusion is based on the con- 
struction of an effective spin-orbit Hamiltonian for the impu- 
rity system of III-V semiconductors, together with its analyti- 
cal and numerical solution. The resulting spin relaxation times 
are in good agreement with the existing experimental values 
for GaAs. The detailed temperature-dependent measurements 
of Ref. il2ll yielded a saturation of Tj below lOK, indicating 
that inelastic processes are irrelevant at low temperatures. We 
thus work with a zero-temperature formalism. 

The envelope-function approximation (EFA) for describing 
conduction-band electrons in zincblende semiconductors in- 



corporates the lattice-scale physics (described by the periodic 
part of the Bloch wave function) into the effective one-body 
Hamiltonian figjlloll 



H - Ho + HsiA + Hbia 
^0 = + V(r) 



2 m* 

i/siA = o- ■ vy X k 



^^BiA = 7 [c^xkAky - k,) + cyclic permutations] . 
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The electrostatic potential V{r) includes all potentials aside 
from the crystal one, while cr is the vector of Pauli matri- 
ces and k = p//i. The effective spin-orbit coupling A, en- 
abled by the structural inversion asymmetry (SIA) is usu- 
ally orders of magnitude larger than the one of vacuum Ao = 
fp-/4tnlc^ ^ 3.7 X 10"'' A^. It can be calculated at the level 
of the 8-band Kane model, which, for example, for GaAs 
yields A ^ -5.3 A^ [20]. The bulk inversion asymmetry 
(BIA) coupling constant y is another material-dependent pa- 
rameter. The exact value of y is a matter of current debate 
ifllH^]- A 14-band model is required for the theoretical esti- 
mation of y, leading to y » 27 eVA^ for GaAs USEIl]. More 
refined theoretical calculations yield somewhat lower values 
II22I - I24II . While early experimental values obtained in bulk 
samples agree approximately with the above-quoted value of 
27eVA^ 12511 . recent results inferred from measurements in 
low-dimensional systems are again consistently lower 

2I. 

In order to study the spin relaxation in the impurity band 
near the MIT, we consider the potential V(r) due to the ionized 
impurities, given by 



y(r) = ;^y,(r) = -;^ 



elr - R„ 



(5) 



where e is the dielectric constant of the semiconductor and 
Rp represents the impurity positions. The potential Vp gives 
rise to hydrogenic states centered at the impurity p. In or- 
der to build the basis of electronic states we will only con- 
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sider the ground state 4>p{r) - <p{\r - Rp|), with 0(r) - where 
(1 / y/na^) exp (-r/a), and a the effective Bohr radius. 

The second-quantized form of the Hamihonian that we 
denote "H, has components 



■TYo = ^ {m' o-\Ho\mo-) c\^,^ c„ 

mi^m' ,c 

- Z 



(6) 

mi^m' ,cr 

'Kso = ^ {m"a\Hso\mo-) c^,_ c,„^ , (7) 

where the label SO stands for SIA or BIA. We denote the 1 s 
state 0„,(r) with spin cr = + 1 in the z-direction by |mcr), and 
cLcr (cmcr) IS the Creation (annihilation) operator of a particle 
in that state (cr = -cr). The matrix elements in Eq. ^ contain 
three-center integrals {m'cr\Vp\mo-} with p + m. Due to the 
exponential decay of 0,„(r), one usually keeps only the term 

R„ 



{m' cr\V,„'\mcr) - -V{ 



J a 



(8) 



with Vo - e^/sa (twice the binding energy of an isolated im- 
purity) and /?,„',„ the distance between the two impurities. The 
resulting Hamiltonian 'Hq defines the well-known Matsubara- 
Toyozawa model (MT) ||29tl . which has been thoroughly stud- 
ied in the context of the MIT. The subtleties, drawbacks and 
applicability of this model to describe the metallic side of 
the MIT, as well as its extension to include the "Hsia spin- 
orbit coupling, have recently been discussed [|30|]. Electron- 
electron interactions induce significant many-body effects on 
the insulating side of the transition, but not on the metallic 
side. Therefore we do not need to include them in our model. 
According to the Mott criterion the critical dimensionless im- 
purity density for the MIT corresponds to TVi = riia^ ^ 0.017, 
corresponding to a critical density of 2 x 10'^ cm"^ for GaAs. 



The matrix element of 'Ksia is 
a-A 

~^ J "' ' |r-R, 

X [(Z - Zm)(rcr - Rm'o-) ' iz - Zm')(fcr ' Rma 



{m' (T\HsiA\mo-) 



dr V(r) ■ 



4>m'{r)4>m{r) 



)], 



(9) 



where rg- - x + icry and R,„,j- - X,n + icrY,n. The Hamilto- 
nian "KsiA represents the generalization of the Rashba cou- 
pling to the case of impurity potentials, and was introduced in 
Ref. H. There, an alternative path to the EFA was followed 
in order to calculate the matrix elements (m'o^|//siAl'wcr), 
which made use of impurity states with spin admixture ob- 
tained from spin-admixed conduction-band Bloch states de- 
rived at the level of the 8-band Kane model. We remark that 
the terms corresponding to p - m, m' in V'(r) give vanish- 
ing contributions to the SIA matrix element due to the axial 
symmetry of the two-center integrals. As a consequence, the 
Rashba matrix elements are given by three-center integrals, 
resulting in very slow spin relaxation OlSll in comparison with 
experimental results. We therefore turn to the Dresselhaus 
term, whose matrix element is given by 

{m' &\HBiA\mo-} = y [{m'\kx(ky - kl)\m) 

+io-{m'\ky(kl-kl)\m}] 

~ T ^y,m'm ~^ ^ ^x.m'ni) ? (10) 



dr- 



e 



-\r-R„ 



a-' J |r - R, 
xia + r)(x - X,„,,„){y 



■m'm\ ^ 
2 



z\ 



(11) 



and L 

m'm is obtained from Ix^m'm with the exchanges Xf^'f^ <— > 
Y,„',n and x «-» y. Performing a rotation of the coordinate sys- 
tem and switching to cylindrical coordinates allows us to eas- 
ily do the angular integral, yielding 



where 



7TC r 

m'm - I dpdzp 

{a+ ^1 



Vp^T(F 



7;o)2(p2+z2)3/2 



^(3z-7?o)-z'(z--Ro) 



c = 2 COS sin Q\\ - sin 0(1-1- sin i^)]. 



(12) 



(13) 



and (7?o, 0, 1^) are the polar coordinates of R^'m in the origi- 
nal reference frame. The integrals in (IT2t are not analytically 
solvable, but they can easily be integrated numerically. This is 
the route that we take later, where we simulate the dynamical 
evolution of initially prepared pure spin states. Before tack- 
ling the numerical problem, we provide a simple estimation of 
the spin lifetime based on the steepest-descent approximation 
of the integral ( fT2t (valid in the limit /?o ^ «)^ given by 

W,«-0«o(a+y)^-*°^^ (14) 

The spatial diffusion of electrons through the network of im- 
purities is accompanied by a spin diffusion characterized by a 
typical spin rotation angle per hop 118 1 



15 Z,„^m' \{m' (T\HBip\m&)'^ 
2 



(15) 



The spin-relaxation time can be defined as the time in which 
the accumulated rotation reaches a unit angle, and therefore is 
given by 

2<02> 
3 ' 



1 



(16) 



where is the typical hopping time. It can be shown that 
this expression is independent of whether the initial state is 
localized or extended y_8]. Taking as the time needed for 
the initial-state population on an impurity site to drop from 1 
to 1/2, we obtain 



1 

Tc 



V2 



^ KmV|//o|mcr)p 



-TV. 



1/2 



(17) 



For the second equality fs?] we have used the impurity aver- 
age assuming a random distribution without hard-core repul- 
sive effects on the scale of the effective Bohr radius lll8ll3lll . 
Performing also the impurity average in Eq. (|T5) we have 

0.258 l-T^) , (18) 
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FIG. 1. Scaling of the spin relaxation time extrapolated to infinite 
system size, with the spin-orbit enhancement factor for densi- 
ties M = 0.02 (circles), 0.029 (squares), 0.037 (triangles), and 0.06 
(crosses). The lines are fits of a quadratic dependence of the relax- 
ation rate on The times are given in units of Ti/Vq. Lower Inset: 
Spin survival probability P for an initial MT eigenstate evolving un- 
der the enhanced Dresselhaus couplings for a density N, = 0.029 
and a system size A' = 3375. Lines of increasing thickness are for 
enhancement factors rj = 75, 100, and 150. Upper Inset: Size depen- 
dence of T^' for Tj = 150. Lines are linear fits to the data that allow 
to extrapolate to the infinite-size values. 



and 



1 



a^Voh ' 



(19) 



This expression for the spin relaxation time is our main re- 
sult. As we discuss below, its confrontation with the exper- 
imental data allows us to identify the Dresselhaus coupling 
as the dominant spin relaxation mechanism in the regime of 
study. Furthermore, numerical calculations of the relaxation 
time within our model provide a complementary path validat- 
ing the analytical approach, since some of the previously used 
approximations can be avoided. 

The numerical procedure starts from the numerical integra- 
tion of ( fTSI ) for a given impurity configuration in order to 
obtain the matrix elements (fTOt for //bia (and similarly for 
Hsia)- We then diagonalize the total Hamiltonian includ- 
ing the two contributions to the spin-orbit coupling, which al- 
lows us to obtain the quantum evolution of an arbitrary state. 
Choosing an initial state with a well-defined spin projection 
(for instance a MT eigenstate with cr = 1) we can follow the 
spin evolution and extract the spin lifetime from it. The weak- 
ness of the spin-orbit coupling translates into spin-admixture 
perturbation energies which are, even for the largest system 
sizes (in terms of number of impurities, A^) that we are able to 
treat numerically, orders of magnitude smaller than the typi- 
cal MT level spacing. This large difference between the two 
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FIG. 2. Spin relaxation time in n-doped GaAs as a function of 
the doping density. The prediction of Eq. ( |19t (solid line) and our 
numerical results (filled circles) for the regime between the metal- 
insulator transition (dark gray line) and the hybridization of the im- 
purity band with the conduction band (light gray line) obtained using 
7 = 27 eVA^ are compared to experimentally obtained values (open 
symbols) using different methods. Circles and squares are from Ref. 
llOll for T = 2K and T = 4.2 K, respectively, along with data from 
Ref. d (diamonds), |7] (triangles), and |12] (reversed triangles). 
The dashed line represents the result of Eq. il9\ with y = 15 eVA^. 



energy scales in finite size simulations masks the spin-orbit- 
driven physics, and forces us to follow an indirect path: we 
introduce an artificially enhanced coupling constant 7/7 and a 
finite-size scaling procedure. The limits — > 00 and then 
77 — > 1 taken at the end of the calculation provide the sought 
spin-relaxation rate. The numerically extracted values of the 
spin relaxation times associated with i/siA are, consistently 
with the results of Ref. lIlSll . considerably larger than the ones 
experimentally observed. Therefore, we hereafter neglect the 
Rashba term in the numerical calculations, and concentrate on 
the spin evolution governed by //bia- 

In the lower inset of Fig. [T]we show typical spin evolutions 
starting from an eigenstate of the MT system with cr = 1 in 
the energy range of extended states of the impurity band, for a 
density TVi - 0.029 just above the MIT transition for three 
values of the coupling constant and - 3375 impurities. 
The initial perturbative regime with a quadratic time decay 
of the spin survival is followed by an exponential decay from 
which the relaxation rate rr' can be inferred, until the sat- 



uration value of 1/2. For each density and effective coupling 
constant rjy the asymptotic value of rj' can be obtained by ex- 
trapolating the finite-A^ values to the infinite size limit (upper 
inset of Fig. [T]). We ran a sufficiently large number of impu- 
rity configurations (typically 40) to make the statistical eiTor 
negligible [33]. The resulting error bars arise from the qual- 
ity of the fittings to the exponential decay of the spin survival. 
In agreement with our analytical results, an inverse quadratic 
dependence of Ts on the coupling strength is obtained (Fig. 
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[T]). Fitting this dependence of Ts on 77 allows us to extract the 
physical values (rj - 1) of Tj. 

In Fig.|2]we present the spin relaxation times resulting from 
our numerical approach for GaAs at four different impurity 
densities above the MIT (black dots), together with the pre- 
diction of Eq. (fTgll (solid line), and the available experimental 
data from RefsHllEI 1211 . The agreement between the ana- 
lytical and numerical approaches is very important in view of 
the complementarity of these two very different ways to ex- 
tract Tj. Both approaches describe the data within the experi- 
mental uncertainty and correctly reproduce the density depen- 
dence of the spin relaxation time. The small departure of the 
analytical and numerical results is not significant, taking into 
account the different approximations of both paths and the ar- 
bitrariness associated with the definition of relaxation times, 
i.e. numerical prefactors in Eqs. (IT6^ and (ITTt . 

While in the critical region and deep into the localized 
regime there is some dispersion of the experimental values 
for GaAs, depending on the different samples and measure- 
ment technique, on the metallic side of the MIT values of Tj > 

1/2 

100 ns are consistently obtained. A decrease of Tj with 
is observed, with a clear change in the density-dependence 
once the impurity and the conduction bands hybridize. Our 
analytical and numerical results of Fig.|2](solid line and filled 
symbols) are obtained using the values Vo = 11 .76 meV and 
y - 27 eVA^ without any adjustable parameter We remark 
that these results are very sensitive to the value of 7. For in- 
stance, taking y = 15 eVA-', that is proposed by some mea- 
surements yl |26|], results in the dashed line. The identifi- 
cation of the DresseUiaus coupling as the dominant channel 
for spin relaxation close to the MIT of 3-dimensional samples 
will help in the precise determination of the material constant 
y, which is crucial for the operation of low-dimensional spin- 
tronic devices. 

In conclusion, we have identified a spin relaxation mech- 
anism characteristic of electrons on the metallic side of the 
metal-insulator transition in the impurity band of semiconduc- 
tors, where up to now a theoretical understanding was lacking, 
thereby solving a longstanding problem in spintronics. Our 
mechanism is derived from the Dresselhaus spin-orbit cou- 
pling. It dominates over the usually stronger Rashba coupling 
in the landscape of hydrogenic impurities in semiconductors 
with zincblende structure, and provides relaxation times that 
agree with the experimentally measured values. 
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